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In 4D compact U(l) lattice gauge theory with a monopole term added to the Wilson action we first reveal some 
properties of a third phase region at negative (3. Then at some larger values of the monopole coupling A by a 
finite-size analysis we find values of the critical exponent v close to, however, different from the Gaussian value. 


1. INTRODUCTION 

In recent higher-statistics studies of 4D com¬ 
pact U(l) lattice gauge theory with the Wilson 
action |jj and with this action supplemented by a 
double charge term || for increasing lattice size 
the critical exponent v has turned out to decrease 
towards i, the value characteristic of a first-order 
transition. Also stabilization of the latent heat 
has been observed. Thus, there are now rather 
strong indications that in those cases the transi¬ 
tion is of first order. In contrast to this, for the 
action where a monopole term with coupling A is 
added, at A = 0.9 the critical exponent has been 
found to be characteristic of second order ||. 

In H the location (3 CI of the transition from 
the confinement phase to the Coulomb phase has 
been determined as a function of A. It has been 
found that this transition line continues to nega¬ 
tive f3. On the other hand, in Refs. |^],[| a further 
transition has been seen at (3 = — 1 for A = 0 and 
at (3 = —0.7 for A = oo. Thus the question arises 
what happens at negative (3 which is addressed in 
Sect. 2. 

The energy distribution indicates that the sec¬ 
ond order region of the confinement-Colomb tran¬ 
sition starts at some finite A above A = 0.7 ||. 
An important question is wether one has a region 
with universal critical properties there. To in¬ 
vestigate this we have performed higher-statistics 
simulations also at A = 1.1 and at A = 0.8 and 
have evaluated the data for several variables by 
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finite-size analyses. The respective results are 
presented in Sect. 3. 

2. PHASE REGIONS 

In |]] the value f3 cr at the transition from the 
confinement phase to the Coulomb phase has 
turned out to decrease with A and to get nega¬ 
tive below A = 1.2 . The symmetry (3 —► — (3, 
Ua —> —Ua of the Wilson action at A = 0 gives 
rise to a transition at f3 = — 1 in addition to the 
one at (3 = 1 |Q. For A y 0 the indicated symme¬ 
try is violated by the monopole term. At A = oo 
only the transition at negative (3 persists and oc¬ 
curs at about [3 = —0.7 [||j|. 

Here we have checked the occurrence of such 
transition at negative (3 also at intermediate val¬ 
ues of A determining the maximum of the specific 
heat C max . It turns out that there is a transition 
line extending from (A, (3) = (0, —1) to (oo, —0.7). 
Since from the properties we have observed so far 
we have no indication of a further subdivision of 
the region below this line, we consider it a third 
phase. Figure 1 gives an overview of the phase re¬ 
gions as they are according to our present know¬ 
ledge. It includes the line separating confinement 
and Coulomb phases obtained in |J and the data 
at negative (3 found here. 

P ne t, the probability to find an infinite network 
of monopole currents (where “infinite” on finite 
lattices is to be defined in accordance with the 
boundary conditions), has turned out to provide 
an unambiguous characterization of the confine¬ 
ment phase and the Coulomb phase |G]. For the 
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Figure 1. Location of phase transition points 0„ 
on 8 4 lattice as function of A between confinement 
and Coulomb phases (circles from C m ax, square 
from P ne t) and to third phase (diamonds). 

periodic boundary conditions used here “infinite” 
means topologically nontrivial in all directions. 
Taking the values 1 and 0 in the confinement and 
Coulomb phases, respectively, P ne t is very effi¬ 
cient to discriminate between those phases. In 
contrast to this in the third phase at fixed (A, 0) 
we have observed values 0 as well as 1 for P ne t- 
That P ne t provides no longer a criterion in the 
third phase can be understood by noting that the 
monopole quantities are not invariant under the 
transformation (3 — » — (3, Ua —■► —Ua- 

As a characteristic feature of the third phase 
we have found that different states exist between 
which transitions in the simulations are strongly 
suppressed. We have observed this phenomenon 
at various negative (3 in the A range from 0 to 2.5 . 
Typical examples of time histories of the average 
plaquette e are given in Figure 2. 

The values of e increase as (3 gets more nega¬ 
tive. They are somewhat below 2 which indicates 
that at sufficiently negative 0 the average value 
of cos 0 gets close to — 1. In view of the the sym¬ 
metry 0 —> —0, Ua —*• —Ua of the Wilson action 
the correspondence of positive cos 0 to positive 0 
and of negative cos 0 to negative 0 is conceivable. 
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Figure 2. Typical time histories of the average 
plaquette e at negative 0 obtained for seven differ¬ 
ent simulation runs, shown for (X, 0) = (2,-20) 
and 8 4 lattice. 

The origin of the different states is not yet clear. 
Such states are similarly observed in spin glasses 
and frustrated systems, and also with sponta¬ 
neous breaking of Z(N) in finite-temperature 
SU(N) gauge theory. 

3. CRITICAL PROPERTIES 

At A = 1.1 and A = 0.8 for each lattize size 
considered Monte Carlo simulations have been 
performed at a number of 0 values in the criti¬ 
cal region. Multihistogram techniques have been 
applied to evaluate the data and the errors have 
been estimated by Jackknife methods. In the 
finite-size analysis in addition to the specific heat 
and the Challa-Landau-Binder (CLB) cumulant 
complex zeros of the partition function, in partic¬ 
ular the Fisher zero zo closest to the 0 axis, have 
been used. 
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Figure 3. Maximum of specific heat C max as 
function of lattice size L for A = 0.8, 0.9 and 1.1 
at transition point /3 cr between confinement and 
Coulomb phases. 


For d = 4 the maximum of the specific heat 
is expected to behave as C max ~ L 3 4 5 6 if the phase 
transition is of first order and as C max ~ L~ if it 
is of second order, where a is the critical exponent 
of the specific heat and v the critical exponent of 
the correlation length. 

In Figure 3 we present the results for C max 
obtained on various lattices. They include data 
from simulations of the present investigation with 
A = 1.1 and A = 0.8 and ones from simulations of 
Ref. U with A = 0.9 . 

The fits to the data in Figure 3 give values for 
^ far from 4, i.e. far from what would be ex¬ 
pected for first order. Using these values and the 
hyperscaling relation a = 2 — d v the values for v 
listed in Table 1 are obtained. They are seen to 
be close to Thus in any case to conclude on 
second order appears quite safe. 

Similar results are obtained for the minimum 
of the CLB cumulant 

^clb = (3T) 

and for the imaginary part of the closest Fisher 


zero Zo . For these quantities finite-size scaling 
predicts the behaviors 

Im(z 0 ) ~ L ~# , (3.2) 

UcLB.min ~ L~~ A . (3.3) 

The results of the respective fits are also listed in 
Table 1. 

The values of v obtained are seen to be close to 
the Gaussian value f, however, different from it. 
The observed increase of v with A could indicate 
a nonuniversal behavior. Another possibility is 
that it is related to finite-size effects. Then the 
increase should disappear on much larger lattices. 
In that case the universal value of v taken on the 
infinite lattice could even be the Gaussian one. 


Table 1 


Critical exponents v from Im(zo), Cy, Uclb- 


A 

Im(zo) 

Cy 

UCLB 

0.8 

0.404(5) 

0.433(2) 

0.421(3) 

0.9 


0.446(5) 


1.1 

0.421(8) 

0.467(2) 

0.455(2) 
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